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ORCA is a low cost remotely operated vehicle which was indigenously 
developed for underwater inspection and survey. As the underwater 
environment is quite unpredictable, dynamic modeling and simulation of the 
remotely operated vehicle are essential to understand the behavior of the 
vehicle and accomplish stabilized navigation. This paper discusses a detailed 
approach to the mathematical modeling of ORCA based on Newtonian 
dynamics and simulating the position and velocity responses in Simulink. 
The open loop nonlinear model of the remotely operated vehicle was used to 
study the navigation challenges due to the various perturbations present 
underwater namely Coriolis and centripetal force, added mass, hydrodynamic 
damping force, and restoring forces. The six-thruster open loop ORCA model 
was subjected to various thrust inputs (25%, 50%, and 75%) to achieve six 
degrees of freedom (DoF) respectively and it was observed that there was 
significant instability in the other DOFs along with the principal direction of 
motion. Further, the authors will incorporate the various control systems in 


ORCA and analyze the stability in navigation induced due to each of them. 
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1. INTRODUCTION 

The tethered underwater robots known as remotely operated vehicles (ROV) are frequently 
employed for underwater survey and inspection. They are unmanned, highly maneuverable underwater 
vehicles that can be used to explore ocean depths while being operated remotely by a user at the water's 
surface. They are extensively utilized for both scientific reasons like ocean exploration and industrial ones 
like internal and external inspections of underwater pipelines and the structural assessment of offshore 
structures. Recently, a number of ROV applications for ocean research have been reported [1]—[4]. Based on 
the application, they range in their size from a micro or mini ROV which is used to reach areas at shallow 
depths to a work class ROV which is used for ocean floor exploration and inspection at unreachable depths. 
They serve to be greatly beneficial due to quick deployment, minimal maintenance, the ability to fit into 
confined areas, extended dive times, and video recording capabilities. Aquaculture, search and rescue, 
underwater intervention, repair, and maintenance work in offshore businesses, such as the oil and gas 
industry, marine structures, marine sciences, naval defense, marine renewable energy, and research reasons 
are just a few of the uses for these underwater ROVs [5]}-[11]. ROVs enables users to capture photos and 
videos to inspect and monitor ports, harbors, and vessels, locate underwater targets, and explore the depths of 
our oceans. ROVs are equipped with manipulators and multiple sensors to conduct survey and inspection 
activities. Sensors like altimeters, pressure sensors, inclinometers, and magnetic and gyrocompass are also 
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used to help ROVs understand their environment [12]. ROVs frequently have cameras for surveillance and 
streaming video. The design and implementation of Aqua bot is discussed [13], an ROV constructed with the 
goal of addressing issues with real-time video streaming and other development methods. Numerous 
localization and navigation techniques, including dead reckoning and inertial navigation systems, are in use. 
These include acoustic navigation, geophysical navigation, fusion of long baseline (LBL), and doppler-based 
techniques [14], [15]. 

In today’s world, computer-aided designs of engineering systems as models are performed. They are 
the next step after designing. They come in handy in documentation and quality assurance. These models can 
hence be used in various ways by firms and people who install, commission, and operate subsea installations 
and equipment. Simulation software has always been a boon to the product development process because 
they give an opportunity to rehearse certain conditions of operation to test for worst-case and best-case 
scenarios. Workflow analyses, contingency plans, and training scenarios for the pilots and technicians who 
will perform the tasks can be thus provided for ease of development. Mathematical models have been 
developed to design the dynamics of an ROV and other underwater vehicles. These models have been 
successfully adapted by many ROVs for performing various activities [16]-[18]. Since an underwater 
channel is highly dynamic, it is necessary to have an effective control system for stabilized movement [19]- 
[21]. Different control systems with varying levels of efficiency can be created for ROVs [22], [23]. 

This paper addresses broadly the development of a mathematical model of the ORCA ROV and 
analyze the position and velocity responses obtained from its simulation. The first section briefs on the 
design and development of ORCA and throws light on its capabilities to conduct an efficient inspection. The 
second section presents detailed mathematical modeling and the necessary design parameters of ORCA 
required to simulate the nonlinear open loop model and study its response. It further projects the Simulink 
model of the ORCA ROV and briefly discusses the various subsystems. Finally, the last section presents the 
performance evaluation of nonlinear open-loop ORCA when subjected to different input thrusts for the 
desired time to study the instability in the position and velocity responses for the six degrees of freedom 
namely surge, sway, heave, roll, yaw, and pitch. 


2. SYSTEM DESIGN 

In this section, the basic design of the ORCA ROV is elaborated. ORCA ROV, as shown in 
Figures 1 and 2, is a 500x350x210 mm? vehicle weighing 10.247 kgs which is equipped with sensors like 
temperature, and pressure, an IMU, and a camera for use in the primary underwater inspection and survey 
application until a depth of 300 m. It was indigenously developed by the authors and team and the novel 
feature of ORCA is shown in Table 1. 


Isometric view Top view 


Figure 2. Fabricated ORCA ROV 
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High-density polyethylene (HDPE) was used to make the six-thruster ROV. It has two body- 
mounted sealed acrylic enclosures for the battery and electronics, respectively. The embedded electronic 
system integrated with sensors as shown in Figure 3 is responsible for achieving the desired movement of the 
vehicle as per the joystick command. The robust navigation is achieved due to the presence of the attitude and 
heading reference system, a sensor fusion algorithm that obtains the IMU raw data to compute the current 
position and attitude of the vehicle. The whole system is powered using the 14.8 V, 18 Ah Li-ion battery. The 
versatility of ORCA has increased due to the freedom in placing the T200 thrusters at two different 
orientations. It has four thrusters mounted at 45 degrees to the x and y axes for active control of surge, sway, 
yaw, and two dedicated thrusters to achieve heave or roll, thereby enabling active control of 5 degrees of 
freedom (DoF) and passive pitch control. An effective thrust vectorization algorithm was developed to 
achieve best-in-class maneuverability. 


Table 1. Novel features of ORCA 


S. No. Unique Features of ORCA Remarks 

1 Flexible thruster mounting. Gives rise to different possible vectored thrust configurations. 
The four horizontal thrusters can be manually swiveled This helps us to customize ORCA as per the underwater 
and mounted at 30, 45 or 60 w.r.t the heading axis. environment to tackle the varying currents. 

2 Open-source software platform (MOOS-IVP) This gives rise to use integrate many sensors to the vehicle in 
Exclusively customized application was developed for ease as the development does not involve any dedicated 
remote operation of the vehicle using this platform. hardware. 

Easy integration of various sensors. 

3 Effective thrust vectorization algorithm Easy and best in class maneuverability with active control in 5 

Algorithm developed on the controller end to achieve degrees of freedom with passive pitch control. 


navigation in 6-DoF. 
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Figure 3. Electronic circuitry scheme in ORCA 


3. OPEN LOOP MODELING OF ORCA 

To understand the motion of ORCA in underwater environments subjected to various perturbations 
it is required to mathematically model the ROV and simulate the position responses. The dynamic model for 
ORCA is derived based on the Newtonian approach as adopted from [21]. The two frames of reference, the 
inertial frame or the earth-fixed frame, and the local body fixed frame with its origin coinciding with the 
vehicle’s center of gravity are considered to define the movement of ORCA ROV in three-dimensional space. 
The assumptions considered to model ORCA ROV involves that i) ROV is a rigid body fully submerged in 
water; ii) disturbance due to a wave is neglected as the ROV is fully submerged in water; and iii) the earth- 
fixed frame of reference is inertial; and ROV is assumed to be moving slowly to carry out the inspection. 

The 6 DoF nonlinear equation of motion of ROV in body fixed frame required to mathematically 
model ORCA is given by (1). 
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Mý + C(v)v +D(v)v+ Ge) = T (1) 


where M is the mass inertia matrix which is the sum of the rigid body inertia mass and added fluid inertia 
mass matrix, M= Mrs +Ma E R“; C(y) is the Coriolis and centripetal force matrix which contains the rigid 
body and added mass terms respectively, C(v) = Cra(v)+Ca(v) E R; D(v) E R®Ć is the damping matrix 
due to the surrounding fluid; Grm) E€ R™! is used to describe the gravitational and buoyancy force on the 
ROV; t € R®! is the total force or thrust vector in a given degree of freedom; v = [ u v w p q r ]" is the 
body-fixed velocity vector and y = [nı n2]" is the earth-fixed position and orientation vector where m1 = [x y 
z ]" is the position vector and n2 = [9 8 y |" is the orientation vector of Euler angles as per Table 2. 


Table 2. Notations used for ROV 


DoF Motion Description Forces and moments ied Se Geer emer ery 
1 Surge (Motion in the x-direction) X u x 
2 Sway (Motion in the y-direction) Y v y 
3 Heave (Motion in the z-direction) Z w Z 
4 Roll (Rotation about x-axis) K p ọ 
5 Pitch (Rotation about y-axis) M q 0 
6 Yaw (Rotation about z-axis) N r y 


Hence, it is evident from (1) that the three external forces, the damping force, the Coriolis force, and 
the gravitational and buoyancy restoring force acting on the ROV influence its navigation in an underwater 
environment. The rigid body mass inertia matrix represents the mass of the ROV and the spread of the mass 
in different directions in terms of the moment of inertia, which is given by 


m 0 0 0 MZg -Myg 
0 m 0 —MZg 0 MX¢ 
0 0 m myg —MXc¢ 0 
Mrs = 0 —MZq myg Ix rly ly 
MZG 0 —MX&G —lyx ly —lyz 
—MYg MXG 0 lex ~ly I, 
where m represents the mass of the ROV, re = [xe ye zeļ! is the center of gravity, and 
I, ly ly 
I= |-lyx ly —lyz| represents the rigid body ROV’s the inertia matrix where I, , Iy and I, are the 
lzx —Izy I, 


moments of inertia about X,Y, and Z axes and the cross products of inertia Ixy = Iyx , Ixz = Lx , and Iyz = Izy. 
ORCA has its center of gravity at re = [xc yc zc] = [0 0 0]" and an inertia matrix of 


0.16419 —0.01493 0.00113 
I= |—0.01496 0.34512 —0.00063 
0.00113 —0.00063 0.24282 


Therefore, the rigid body mass matrix of ORCA is 


10.247 0 0 0 0 0 
0 10.247 0 0 0 0 
TPE 0 0 0 0 0 0 
RB 0 0 0 0.16419 —0.01493 0.00113 
0 0 0 —0.01493 0.34512 —0.00063 
0 0 0 0.00113 -—0.00063 0.24282 


The inertial force known as the Coriolis force acts on an object moving within a frame of reference that 
rotates with regard to another object. The Coriolis and centripetal forces on ROV are given by matrix as per 
the Newtonian dynamics: 


m [vz XV, + V2X(V2XIcg)] m [vz xvi + V2 xX (V2 XI) | 


Cre = v2 X (1v2) + mrg X (v2 X V1) v2 X (Iv2) + mrg x (V2 X v1) 
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where m is the mass of ROV, vı = [u v w]" and v2=[p q r]” are the velocity vectors, re is the center of gravity 
of the ROV and I is the moment of inertia of the vehicle. The motion of the fluid around the vehicle due to 
the motion of the ROV causes a hydrodynamic negating force which will retard the motion of the ROV. The 
added mass matrix is given by: 


Xa 0 0 0 0 0 

oY% 0 0 0 0 

0.0 Zy > “0° 0 
Ma=!0 0 0 Ka O 0 
0 0 0 0 My 0 

0 0 0 0 0 % 


where Xa Yy, Zw, Ky, Ma and N; are the hydrodynamic derivatives in the respective DoFs. The added mass 
matrix of ORCA is given by: 


Ma = —diag{ 2.98509, 7.30064, 17.23463, 0.38282, 1.87219, 0.719066} 


The corresponding Coriolis and centripetal added mass matrix Ca(v) are given by: 


0 0 0 0 -Zw Yv 
0 0 0 Zw 0 -Xw 
0 0 0 -Yv Xu 0 
C=] 0 -Zw Yv 0 -Nr Maq 
Zw 0 -Xu Nr 0  —Kpp 
-Yv Xau 0  -Maq Kọp 0 
0 0 0 0 ~17.23463w 7.30064v 
0 0 0 17.23463w 0 ~2.98509u 
0 0 0 ~7.30064v 2.98509u 0 
= 0 ~17.23463w 7.30064v 0 -0.719006r  1.87219r 
17.23463w 0 —2.98509u 0.719006r 0 —0.38282p 
~7.30064v —2.98509u 0 ~1.87219r  0.38282p 0 


Hydrodynamic damping appears as the sum of the linear damping caused due to skin friction and quadratic 
damping due to vortex shedding. Hence, D(v) = Dy, + Dg, where 


Dy = —diag{ Xu, Yv, Zw, Kp, Mg, Ne} 
De = - diag {Xujup, Yvivp Zwlwp Kpip Malar Nriris 


represents the linear and quadratic damping coefficients. Calculating the respective drag coefficients for 
ORCA in 6-DoF, we obtain 


D; = —diag{0.1682, 0.3727, 2.7871, 0, 0, 0} 
Do = —diag{ 11.612, 46.197, 86.647, 0, 0, 0.3168} 


The restoring gravitational and buoyancy force vector acting on the ROV is given by 


(W — B) sin O 
—(W — B) cos 8 sin ọ 
Gf (n) = —(W — B) cos 8 cos b 
—(ycW — ygB) cos 8 cos p + (ZW — ZgB) cos 8 sind 
(ZGW — ZgB) sin 8 + (xW — xgB) cos 8 cos p 
—(xcW — xgB) cos 8 sin p — (yW — ygB) sin 0 


where W is the weight of the ROV, B represents the ROV buoyancy, and [¢ 8 J] represents the Euler angles 
representing the orientation of the ROV. ORCA is designed to be positively buoyant with center of 
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Buoyancy, rg = [Xg Yg Zg]” = [ 0 0 -0.0082] T. The weight and buoyancy of ORCA are found to be 
100.52307 N. Hence the buoyancy and gravitational force matrix are 


0 
0 
0 
0.8242 cos 8 sind 
0.8242 sin 0 
0 


Gf (n) = 


The input thrust generated by each of the six thrusters, u = [u4 Uz U3 U4 Us Us] is related to the 
force and moment vector generated in the 6-DoF, T = [Tx Ty Tz Tp To Tyl” as 


t= Tu (2) 


where T is the thrust configuration matrix which represents the position of each thruster from the center of 
gravity (COG) of the vehicle. Four of the six thrusters of ORCA are mounted at a 45° to the x and y-axes 
about the z axis, and the other two thrusters are mounted with its axis parallel to the z-axis. The thrust 
vectorization of ORCA is shown in Figure 4. The force and moment vector generated in 6-DoF is given by 
(3) to (8). 


u; sina + uz sina + u,sina + ussina = Ty (3) 
U4 COS Q — Uz COS Q — U4 COS Q + Us COSA = Ty (4) 
uz + U4 = T, (5) 
—yu; + yu, = Te (6) 
u,Zsina+ uzsin a + uzZz sin a + u4zZsina = Te (7) 
U1X COS A — U2X COSA + UZXCOSA — U4% COSA = Ty (8) 


Figure 4. Thrust vectorization of ORCA 


Hence, the thrust configuration matrix for ORCA is given by 


sina sina 0 0 sina sina 
cos a —cos a 0 0 -cosa cos a 
z 0 0 1 1 0 0 
T=] o 0 xx oO 0 
zsina zsina 0 0 zsina zsina 
ycosa —ycosa 0 0 ycosa -ycosa 


where @ is angle at which the thruster is mounted (30°, 45°, or 60°), x is distance of the thruster from the 
COG along x axis, y is distance of the thruster from the COG along y axis and z is the distance from the COG 
along z axis. By substituting: 
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Tx 0.7071 0.7071 0 0 0.7071 0.70717 pus 
Ty 0.7071 —0.7071 0 0 -0.7071 0.7071 | |U2 
T| 0 0 1 1 0 0 U3 
To] 0 0 —0.110 0.110 0 0 u, 
To 0.0176 0.0176 0 0 0.0176 0.0176 | |Us 
Ty 0.1237 —0.1237 0 0 0.1237 —0.12371 lug 


the dynamics in the body-fixed frame are transformed to earth-fixed frame using Euler’s angles. The 
kinematic equation representing the Euler transformation is (9), 


n=I@Mv (9) 


where the Euler transformation matrix J(n) is taken from a successive rotation of the Euler’s angles n, = 
[p@w]? about Z, Y, and X axes. 


l= lo sasl K 


1) = PD” panl 


where 


cpc) —sCh)c(h) + cQ)s(8)s(h) — s(p)s(h) + cC)c(H)s(8) 
J12) = |s@)c(8)— capele) + s(H)s(@)s() — —c()s() + s(©)s(p)c(0) 


—s(0) c(8)s(#) c(8)c() 
1 s(h)t(®) c(p)t(8) 
0 c(d) —s() 
Ja2) = sb) c() 
° @ @ 


where s=sin(.), c=cos(.), and t=tan(.). Therefore, this transformation converts the velocity vector in the 
body-fixed frame to the position and orientation vector in earth fixed inertial frame. Having the model as the 
basis, the simulation of the ORCA ROV was performed in Simulink to analyze the position and velocity 
response for the 6-DoF and understand the navigation challenges in the present open loop system. 

As shown in Figure 5, the nonlinear open loop model consists of five subsystems namely, thrust 
configuration block, inertial and added a mass block, Coriolis and centripetal block, gravitational force block, 
and Euler transformation block. The nonlinearity in the model is due to the external forces impacting its 
navigation in the underwater environment. This model with the subsystems realizes (10). 


v= MHri- [C(v)v + D(v)v + Gm) (10) 


The acceleration vector obtained is integrated to form the body-frame velocity vector and then 
further transformed into the inertial frame position and orientation vector. The position and velocity 
responses obtained for a given input thrust and desired simulation time are analyzed in the following section 
to understand the navigation challenges in ORCA ROV. 


4. PERFORMANCE ANALYSIS OF NONLINEAR OPEN LOOP POSITION AND VELOCITY 
RESPONSES FOR ORCA 

The performance analysis involved a lot of intensive testing of the model of ROV-ORCA, by 
varying the input parameters and observing the outputs. ORCA is a six-thruster model with thrusters 1, 2, 5, 
and 6 and thrusters 3 and 4 as shown in Figure 5 for movement in the horizontal plane and the vertical plane, 
respectively. The authors have considered a shallow water channel for modeling which has a density of 
998.2 kg/m?, viscosity of 0.001003 N.s/m?, and pressure of 1.02 bar at 25 °C. Thrust for each of the 6 
thrusters is given in the form of input and the navigation behavior of ORCA is observed at the output by 
means of velocity and position response graphs. ORCA is equipped with T200 thrusters which can produce a 
maximum of 4.53 KgF in a clockwise direction and 3.5 KgF in an anticlockwise direction. Varying input 
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thrust was provided, and the responses were generated with input to be 25% of the maximum thrust advisable 
for the thrusters, then at 50% thrust, and then at 75% thrust. At each of these input cases, a time variation 
from 10 seconds to 300 seconds was given, and the responses were recorded. Shown here are responses to an 
input of 50% thrust (2.265 KgF). 


Output (Position) 
| eout | 


Euler Transformation 


acooaraðon 
pos 


added mass 
Corclts & Centripetal vel 
gravitational force 


MA 


Desired Input(Thrust) 
L 
TI(LF) 
| 
TARF} 


Figure 5. Nonlinear open loop ROV model for simulation 


Modeling and performance analysis of open-loop remotely operated vehicles ORCA (Tejaswini Panati) 


116 o ISSN: 2722-2586 


Figure 6 represents the velocity and position response of ORCA when commanded to surge forward. 
To achieve this response, the thrusters T1, T2, T5, and T6 in Figure 5 were provided with +2.265 KgF 
thrust for 300 seconds. As observed in Figure 6, ORCA actively surged with a velocity of 0.6-0.8 m/sec 
to cover a distance of 200 m in 5 minutes. Surge being the principal direction of navigation in this case 
there is a predominant increase in the distance moved by ORCA in the x direction. While it can be observed 
that in the other linear directions like the sway and heave there is still some amount of disturbance due to 
the action of the damping and Coriolis forces due to which the ROV ends up at 1m to the right and 0.4 m 
deep after 300 seconds. While the disturbance in the form of roll and pitch is minimal it is observed that 
the system yaws to almost 4.5 degrees. Hence, there is some instability induced in other degrees of freedom 
along with the movement of the ROV in the x direction. 
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Figure 6. Performance of ORCA in surge direction: (a) velocity and (b) position response 
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Figure 7 represents the velocity and position response of ORCA when commanded to sway right. To 
achieve this response, as per the thrust vectorization, the thrusters T1 and T6 in Figure 5 were provided with 
+2.265 kgf thrust, and T2 and T5 were provided with -2.265 kgf thrust for 300 seconds. As observed in 
Figure 7, ORCA sways right with a velocity of 0.3-0.4 m/sec to cover a distance of 100 m right in 5 minutes. 
Along with swaying right, it is observed that there is still some amount of disturbance due to the action of the 
damping and Coriolis forces due to which the ROV ends up at 2.2 m forward and 0.4 m deep after 300 


seconds. In comparison to the surge case, the roll induced in the vehicle has been reduced while the pitch and 
yaw effects remain the same. 
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Figure 7. Performance of ORCA in sway direction: (a) velocity and (b) position response 
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Similarly, for each degree of freedom (heave motion in Figure 8, yaw motions in Figure 9, and roll 
motions Figure 10), a careful analysis was performed and keen observations on other non-principal directions 
were made. Hence, an overall inference from the performance analysis is the thrust input given to each 
thruster increases for a given DoF, the distance moved in the principal direction for the same time interval 
also increases. It is observed that with time, the error gets aggregated. Along with the negating forces 
(damping, Coriolis, and restoring force) acting on the ROV, it reflects as a disturbance in other directions too. 
This gave the authors a clear picture of instability in certain degrees of freedom. These directions are 
recommended to use a controller for smooth navigation as it helps us in the reduction of errors. Hence, the 
control system plays a vital role in nullifying unnecessary deviations and helps us achieve stabilized 
navigation. 
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Figure 8. Performance of ORCA in heave direction: (a) velocity and (b) position response 
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Figure 9. Performance of ORCA in yaw direction: (a) velocity and (b) position response 


ORCA ROV post-fabrication and successive bench tests were deployed in a shallow water region, as 
shown in Figure 11, to study the behavior of the open loop system (12.931420858589528, 80.22865283652679). 
Multiple straight-line missions were completed using an open-loop ORCA for a time span of 5-15 mins each 
along with the IMU connected. The data gathered during these missions were represented graphically as 
shown in Figures 12 and 13 for comprehensive analysis of the effect of various environmental parameters on 
the vehicle as well as to study its dynamics as it moves. 
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The data from IMU was recorded for two successful missions or trials. Mission 1, as shown in 
Figure 12(a) and (b), focused on forward surge movement where it does show a lot of deviations from the 
principal direction of movement. This is seen as spikes in the graph. Mission 2, as shown in Figure 13 (a) and 
(b), also has similar spikes indicating sudden and minor deviation of movement of the vehicle from its 


principal direction. 
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Figure 10. Performance of ORCA in roll direction: (a) velocity and (b) position response 
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Figure 11. Real time deployment of ORCA 
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Figure 12. IMU post-processed graphs for (a) the angular velocity versus time and (b) acceleration versus 
time in mission 1 
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Figure 13. IMU post-processed graphs for (a) acceleration versus time and (b) angular velocity versus time in 


mission 2 


5. CONCLUSION 
This paper aims at discussing the effect of an open-loop controller on the position and velocity 
responses of ROV-ORCA. The various response graphs suggest that, given 50% thrust to the respective 
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thrusters to achieve a particular direction of motion, the movement of ROV is not only observed in the 
principal direction but also in the other directions. Although insignificant, this slight deviation from the path 
suggests that as the settling time for ROV motion increases, not only does the distance moved from the initial 
position increase, but also the accumulated error due to the dynamics of the model also increases. 

As shown, the open loop ROV induces a significant amount of instability in the other degrees of 
freedom along with the principal direction of motion. So, the current scope of the work involves analyzing 
and incorporating various control systems which will enhance the navigation of the vehicle underwater and 
reduce the error or deviation induced by the external perturbations present underwater. Further, an intense 
performance analysis is performed on this nonlinear closed-loop ORCA ROV to study its controlled behavior 
in response to the input command. The control system with a good stability score will be implemented in 
ORCA which will receive feedback from the inertial measurement unit about its current state and reduce the 
error to achieve stabilized navigation. 
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